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ABSTRACT 

We  show  uhat,  under  reasonable  assmrptions,  any  collijion- 
avoiding  motion  planning  problem  for  a  moving  system  with 
two  degrees  of  freedom  can  be  solved  in  time  0(\,(n)  log^n), 
where  n  is  the  number  of  collision  constraints  imposed  on  the 
system,  5  is  a  fixed  parameter  depending  e.g.  on  the  maximum 
algebraic  degree  of  these  constraints,  and  X^(n)  is  the  (almost 
linear)  maximum  length  of  {n,s)  Davenport  Schinzel  sequences. 
This  follows  from  an  upper  bound  of  0(X,(/i))  that  we  establish 
for  the  combinatorial  complexity  of  a  single  connected  com- 
ponent of  the  space  of  all  free  placements  of  the  moving  system. 
Although  our  study  is  motivated  by  motion  planning,  it  is  actu- 
ally a  study  of  topological,  combinatorial,  and  algorithmic  issues 
involving  a  single  face  in  an  arrangement  of  curves.  Our  results 
thus  extend  beyond  the  area  of  motion  planning,  and  have  appli- 
cations in  many  other  areas. 


1.  Introduction 

Let  fi  be  a  robot  system  having  two  degrees  of  freedom  (a  2-DOF  system 
in  short),  which  is  free  to  move  in  some  two  or  three  dimensional  space 
amidst  a  finite  set  of  m  pairwise  openly  disjoint  obstacles  Oi,  .  .  .  ,0„, 
whose  geometry  is  known  to  the  system.  The  set  of  all  placements  of  fl  is  a 
2-D  parametric  space  which  we  denote  by  AP  (the  space  of  "all  placements"; 
intuitively  we  think  of  AP  as  a  locally  Euclidean  manifold  which  can  be 
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covered    by    a    constant    and    usually    small    number    of    patches,    each 
homeomorphic  to  some  simple  planar  domain) .  The  subset  of  AP  containing 

m 

all  free  placements  of  B  (i.e.  placements  in  which  B  does  not  intersect  U  Oj, 

and  no  two  subparts  of  B  intersect  one  another)  is  referred  to  as  the  free  con- 
figuration space  of  B,  and  is  denoted  by  FP  (the  space  of  "free  placements"); 
we  also  denote  by  BFP  the  (topological)  boundary  of  FP.  The  motion  plan- 
ning problem  for  B  is:  given  an  initial  placement  Zi  and  a  (desired)  final 
placement  Zi  of  B,  both  belonging  to  FP,  determine  whether  there  exists  a 
collision-free  continuous  motion  of  B  between  these  placements,  i.e.  a  con- 
tinuous path  from  Zi  to  Z2  within  FP  (often  we  are  willing  to  "compromise" 
and  require  only  that  the  path  lie  in  FP  U  BFP).  If  so,  we  also  wish  to  plan 
such  a  motion.  Since  for  such  a  motion  to  exist,  the  two  given  placements 
must  lie  in  the  same  connected  component  of  FP  (or  of  FP  U  BFP),  an 
equivalent  formulation  of  the  problem  is  to  calculate  (a  discrete  combinatorial 
representation  of)  the  decomposition  of  FP  (or  of  its  closure)  into  (arcwise) 
connected  components. 

This  is  a  relatively  simple  instance  of  the  general  motion  planning  prob- 
lem; see  below  for  a  review  of  existing  work  on  the  problem. 

In  our  case,  the  boundary  BFP  of  the  space  FP  can  in  general  be  defined 
by  a  finite  number  n  of  collision  constraints  y-i,  .  .  .  ,"y„,  each  of  which  is  the 
locus  of  placements  of  B  at  which  contact  is  made  between  a  specific  subpart 
of  B  and  a  specific  subpart  of  some  obstacle  or  between  two  specific  subparts 
of  B.  We  will  assume,  as  is  customary  in  motion  planning,  that  these  con- 
straints are  connected  and  simple  algebraic  arcs  of  some  small  and  fixed  max- 
imal degree  d  (but  the  results  obtained  in  this  paper  also  hold  under  weaker 
assumptions  -  see  below). 

With  no  real  loss  of  generality,  we  assume  the  space  AP  of  all  system 
placements  to  be  planar.  In  general,  the  manifold  AP  can  be  more  complex. 
For  example,  for  a  planar  "Stanford  arm",  namely  a  line  segment  free  to  slide 
through  a  fixed  point  and  also  to  rotate  around  it  [FWY],  AP  is  cylindrical; 
for  a  two-link  planar  robot  arm,  AP  is  toroidal,  etc.  In  these  cases  we  break 
AP  into  0(1)  planar  patches,  apply  our  analysis  in  each  patch  separately,  and 
then  glue  the  results  together. 

Let  us  form  the  arrangement  A  of  the  arcs  -yi,  .  .  .  ,•>/„.  This  is  the 
planar  map  whose  vertices  are  the  endpoints  of  these  arcs  and  all  their  inter- 
section points,  whose  edges  are  maximal  connected  subarcs  of  the  -y/'s  whose 
interiors  do  not  meet  any  other  arc,  and  whose  faces  are  the  connected  com- 

n 

ponents  of  the  complement  of    U  7/.    By  Bezout's  Theorem,   the  total 
number  of  vertices  of  A  is  0{d'^n^)=0{n^);  it  easily  follows  from  Euler's 


formula  that  the  number  of  edges  and  faces  of  A  is  also  0{n^). 

The  above  c»nsiderations  imply  that  each  connected  component  of  FP 
must  be  a  face  of  A.  The  converse  is  not  true  in  general,  because  some  faces 
of  A  might  represent  regions  of  forbidden  placements  of  B.  Nevertheless, 
once  A  is  available,  we  can  obtain  FP  from  it  in  a  straightforward  way,  by 
pruning  away  the  forbidden  faces.  In  particular,  it  follows  that  the  combina- 
torial complexity  of  FP  is  0(n^).  Calculation  of  A  can  be  accomplished  by  a 
standard  line  sweeping  technique  (as  in  [BO]),  in  time  0((n+p)  log  n), 
where  p  is  the  total  number  of  intersections  between  the  curves  7/,  and  is 
O(n^)  in  the  worst-case.  An  alternative  technique  has  been  recently  pro- 
posed in  [EGPPSS];  it  runs  in  sUghtly  more  than  quadratic  time,  and  will  be 
mentioned  below.  Another  related  technique  is  given  in  [Q]. 

So  far  we  have  roughly  quadratic-time  methods  for  producing  FP.  More- 
over, these  methods  are  close  to  optimal  in  the  worst-case,  because  there  are 
many  cases  of  motion  planning  problems  with  two  DOFs  in  which  FP  is  actu- 
ally of  quadratic  size  (see  below).  The  key  observation  of  the  present  paper 
is  that  in  most  cases  one  does  not  need  to  calculate  the  entire  FP,  but  only  its 
connected  component  C  containing  the  initial  given  placement  Zq  of  B. 
Indeed,  as  long  as  B  moves  in  a  collision-free  manner  from  Zq  (and  is  not 
artificially  "lifted  up"  and  "re-started"  in  a  position  lying  in  a  different  com- 
ponent of  FP),  it  will  have  to  remain  within  C.  Our  goal  is  thus  to  pre- 
calculate  only  the  component  C,  rather  than  the  entire  FP,  and  thereby 
achieve  better  performance. 

Our  first  main  result  is  that  the  combinatorial  complexity  of  such  a  single 
connected  component  C  is  only  C>(X,+2(n)),  where  s^d^  is  the  maximum 
number  of  intersections  of  any  two  constraint  curves  7/,  yj.  Here  X;.(n) 
denotes  the  maximum  length  of  {n,r)  Davenport-Schinzel  sequences,  i.e. 
sequences  composed  of  n  symbols,  which  do  not  contain  equal  adjacent  ele- 
ments and  also  do  not  contain  an  alternating  subsequence  of  two  distinct  sym- 
bols of  length  r  + 2.  It  is  known  ([HS],  [ASS])  that  k,{n)  is  almost  linear  in  n 
for  any  fixed  r.  More  specifically, 

Xi(n)  =  n;  Xiin)  =  2n-l  (trivial). 

X3(n)  =  0(na(n)),  where  a(n)  is  the  functional  inverse  of  Ackermann's 
function,  and  thus  grows  extremely  slowly  [HS]. 

X4(n)=  e(n-2«(''))[ASS]. 

X2.(n)  =  0(n-2^(«('')'"')),  for  5>2  [ASS]. 

X2,  +  i(n)  =  0(na(n)0(«('')'"')),  for  j>2  [ASS]. 

\2sin)  =  n(«-2"(°('')'"')),  for5>2  [ASS]. 


Our  result  is  actually  a  general  topological  property  of  arrangements  of 
curves.  That  is,  a  single  face  in  such  an  arrangement  of  n  arcs  in  the  plane, 
any  two  of  which  intersect  in  at  most  s  points,  has  combinatorial  complexity 
0{ks+2{^))-  This  result  has  already  been  established  in  [PSS]  for  the  special 
case  of  arrangements  of  line  segments,  where  s=l  and  the  complexity  of  a 
single  face  is  thus  0(X3(n))  =  0(na(n)).  Our  result  extends  a  recent  (some- 
what simpler)  result  in  [SS5],  showing  that  if  the  -y,'s  are  closed  Jordan 
curves  then  the  complexity  of  a  single  face  of  A  is  only  0(X,(n)). 

Our  second  main  result  is  an  algorithm  which,  given  a  collection  of  n 
arcs  7i,  .  .  .  .-yn,  and  a  point  Zq  not  lying  on  any  of  them,  calculates  the  face 
of  the  arrangement  of  these  arcs  which  contains  Zq.  Here  we  assimie  that 
any  two  arcs  -y^,  7^  intersect  in  at  most  s  points,  and  that  any  arc  -y/  has  at 
most  t  points  of  vertical  tangency  (so  that  it  can  be  broken  up  into  at  most 
/  +  1  ;c-monotone  subarcs);  here  s  and  t  are  assumed  to  be  small  fixed  con- 
stants. Moreover,  we  assume  a  model  of  computation  where  certain  primitive 
operations  involving  one  or  two  arcs  take  constant  time;  typical  such  opera- 
tions are:  finding  the  intersection  points  of  a  pair  of  arcs,  finding  the  points 
of  vertical  tangency  of  a  given  arc,  finding  the  intersections  of  an  arc  with  a 
vertical  line,  etc.  All  these  assumptions  are  reasonable  if  the  arcs  y,  are  alge- 
braic of  low  degree  as  assumed  above.  Under  these  assimiptions,  our  algo- 
rithm nms  in  0(Xj+2(n)  log^n)  time  and  0(X,+2('»))  space.  Again,  in  the 
special  case  discussed  in  [PSS]  (where  the  7/'s  are  line  segments  arising  as  the 
Minkowski  difference  of  the  boundaries  of  two  simple  polygons),  another 
algorithm,  of  comparable  complexity,  has  been  presented.  Recently, 
Edelsbrunner,  Guibas  and  Sharir  [EGS]  have  obtained  an  0(na(n)  log^n) 
algorithm  for  calculating  a  singe  face  in  an  arbitrary  arrangement  of  line  seg- 
ments. We  extend  the  technique  of  [EGS]  to  handle  arrangements  of  curved 
arcs.  The  previous  algorithm  in  [PSS]  is  based  on  ray  shooting  in  simple 
polygons  (see  [CG]).  This  technique  can  be  applied  only  in  very  restricted  cir- 
cumstances (some  more  of  which  are  mentioned  below  and  discussed  in  [Si] 
in  more  detail);  in  particular,  it  does  not  seem  to  generalize  to  curved  arcs. 
The  technique  of  [EGS]  is  based  on  line  sweeping,  is  conceptually  simpler, 
and  can  be  extended  to  our  case,  although  this  generalization  requires  some 
additional  analysis  of  the  intersection  pattern  of  such  arcs,  which  is  developed 
below. 

The  special  cases  (in  addition  to  that  in  [PSS])  in  which  the  ray  shooting 
technique  can  be  appHed  include: 

(a)  Translational  motion  of  a  simple  polygon  P  with  k  sides  amidst  a  collec- 
tion of  n  point-obstacles  (imagine  P  translating  on  a  board  amidst  a  col- 
lection of  pins  or  pegs  tacked  to  the  board) .  Here  we  can  calculate  a  con- 
nected component  of  FP  in  time  O  (nik  a  (nit)  log  nit  log  n)   (which  is 


slightly  better  than  our  general  bound). 

(b)  The  case  in  which  the  forbidden  subspaces  induced  by  the  problem  con- 
straints are  all  polygonal,  and  their  sides  have  only  a  fixed  and  small 
number  of  possible  orientations.  This  would  be  the  case,  e.g.  for  transla- 
tional  motion  of  a  rectilinear  simple  polygon  amidst  a  coUection  of  rectil- 
inear obstacles.  In  this  case  we  can  calculate  a  single  connected  com- 
ponent of  FP  in  time  0{nk  a(n)  log  n),  which  reduces  to  0{n  log  n)  in 
case  of  rectilinear  regions,  and  which  again  is  better  than  our  general 
bound. 

However,  as  it  turns  out,  these  time  bounds  can  also  be  obtained  by  an 
appropriate  fine-tuning  of  the  recursion  of  the  general  algorithm  given  in  this 
paper.  For  this  reason  we  omit  here  details  concerning  this  alternative  ray- 
shooting  technique. 

Some  related  results  for  the  case  of  lines  or  line  segments  are  given  in 
[EGHSSSW],  where,  for  example,  it  is  shown  how  to  preprocess  a  collection 
of  n  lines  in  (randomized)  time  roughly  0{n^'^)  and  roughly  linear  space  so 
that  the  face  in  the  arrangement  of  these  lines  containing  any  given  query 
point  can  be  retrieved  in  time  roughly  0{n^'^  +  k),  where  k  is  the  size  of  the 
face. 

As  mentioned  in  the  abstract,  we  regard  our  topological  and  algorithmic 
results  as  basic  important  properties  of  arrangements  of  curves  in  the  plane. 
Our  results  have  recently  been  applied  to  various  other  problems.  One  appli- 
cation in  [EGPPSS]  obtains  a  generalized  "horizon  theorem"  for  arrange- 
ments of  curves,  and  an  incremental  algorithm  for  constructing  such  arrange- 
ments in  roughly  quadratic  time.  Another  recent  application  in  [AS]  obtains 
fast  algorithms  for  detecting  intersections  between  two  collections  of  arcs  in 
the  plane.  Our  results  have  also  been  recently  applied  in  [SSi]  to  obtain  effi- 
cient coordinated  motion  planning  algorithms  for  two  independent  systems 
with  two  degrees  of  freedom  each. 

Related  work 

In  an  initial  series  of  papers  Schwartz  and  Sharir  ([SSI],  [SS2],  [SS3], 
[SA],  [SS4])  obtained  polynomial-time  motion  planning  algorithms  for  the 
general  algebraic  case  and  for  several  specific  robot  systems.  The  Schwartz- 
Sharir  algorithms  involve  decomposition  of  FP  into  finitely  many  simple  con- 
nected cells  and  construction  of  a  connectivity  graph  CG  representing  adja- 
cency of  these  cells  in  FP.  Several  recent  improvements  of  these  initial  results 
have  been  based  on  generalized  Voronoi  diagrams  ([OY],  [OSYl],  [0SY2], 
[LS2]),  whereas  others  involve  optimized  variants  of  the  cell  decomposition 
method  ([LSI]).   Another  recent  technique,  due  to  Sifrony  and  Sharir  [SiS], 


obtains  a  motion  planning  algorithm  for  the  special  case  of  a  line  segment  (a 
"rod")  translating  and  rotating  in  a  2-D  polygonal  region  by  an  explicit  calcu- 
lation of  the  boundary  of  FP. 

A  special  case  of  the  2-DOF  motion  planning  problem  is  that  of  planning 
a  purely  translational  motion  for  a  planar  object  B  amidst  polygonal  obsta- 
cles. This  problem  has  been  studied  in  [OY]  for  the  case  where  fi  is  a  disc, 
and  in  [KS],  [KLPS],  [BZ]  and  [LS2]  for  the  case  where  B  is  a  convex 
polygon.  These  purely  translational  cases  turn  out  to  be  more  favorable  than 
general  2-DOF  problems,  in  that  the  FP  boundary  in  each  of  the  above  two 
cases  contains  only  0{n)  vertices,  and  can  be  optimally  calculated  in 
0(n  log  n)  time  using  generalized  Voronoi  diagrams  (cf.  [OY],  [LS2]).  In 
general,  however,  the  FP  boundary  for  2-DOF  motion  planning  problems  can 
contain  fl(n^)  vertices  (this  happens  even  in  the  purely  translational  case 
when  the  moving  system  is  a  non-convex  polygon;  see  a  remark  in  [KS]). 

Additional  discussion  of  related  results  is  given  in  the  following  section 
2. 

The  paper  is  organized  as  follows.  In  section  2  we  introduce  the  termi- 
nology and  derive  a  few  initial  observations  about  the  problem  structure.  In 
section  3  we  analyze  the  combinatorial  complexity  of  a  single  component  of 
FP,  or,  more  generally,  of  a  single  face  in  an  arrangement  of  curves.  In  sec- 
tion 4  we  present  an  efficient  algorithm  for  calculating  a  single  sudi  face,  and 
in  section  5  we  conclude  with  a  discussion  of  several  simple  cases,  extensions, 
and  open  problems. 

2.  Terminology  and  Initial  Analysis 

As  already  mentioned,  we  adopt  the  following  abstract  representation  for 
the  general  2-DOF  motion  planning  problem.  For  simplicity  of  exposition, 
we  assume  that  the  space  AP  of  all  possible  (not  necessarily  free)  configura- 
tions of  the  moving  system  B  can  be  embedded  in  the  Eudidean  plane.  The 
space  FP  is  assumed  to  be  defined  as  a  Boolean  combination  of  n  algebraic 
inequalities  whose  maximal  degree  is  at  most  some  fixed  (usually  small)  con- 
stant d.  In  the  simple  case  at  hand,  it  is  easily  checked  that  the  boundary  of 
FP  is  contained  in  the  union  of  0(n)  (algebraic)  simple  Jordan  arcs 
"Yi.  •  •  •  .7n»  where  each  -y,  is  the  locus  of  placements  Z  of  B  in  which  a 
specific  subpart  of  B  touches  a  specific  part  of  some  obstacle,  or  two  subparts 
of  B  touch  one  another  (if  this  locus  is  not  a  simple  curve,  we  break  it  into  a 
small  number  of  simple  connected  arcs).  Consider  the  arrangement  A  of  the 
curves  -y/,  i.e.  the  planar  map  whose  vertices  are  the  endpoints  and  intersec- 
tion points  of  the  curves  -y/,  whose  edges  are  the  connected  components  of 
the  curves  yt  when  these  vertices  are  removed,  and  whose  faces  are  the 


connected  components  of  the  complement  of  the  union  of  these  curves.  Thus 
each  face  of  A  either  consists  of  only  free  placements  of  B  (and  is  thus  a  con- 
nected component  of  FP),  or  consists  of  only  "forbidden"  placements  (where 
actual  penetration  of  a  part  of  B  into  an  obstacle,  or  of  one  part  of  B  into 
another,  takes  place).  Let  Zq  iFP  be  a  given  initial  placement  of  B.  As  above, 
our  goal  is  to  calculate  the  connected  component  of  FP  containing  Zq,  i.e.  the 
face  of  A  containing  Zq  . 

When  AP  cannot  be  globally  embedded  in  the  plane,  e.g.  when  one 
degree  of  freedom  6  is  rotational  and  admits  the  full  2ir  range  of  orienta- 
tions, we  will  represent  AP  as  the  union  of  finitely  many  planar  "patches"  - 
in  the  above  example,  these  would  correspond  to  the  subranges  0<e^'ir  and 
'ir<0<2'iT  -  and  apply  the  foregoing  analysis  to  each  of  them  separately.  It 
"s  importan*  to  note  Ihct  the  uumber  of  patches  is  independent  of  the 
geometric  complexity  of  the  workspace  of  the  system  B,  and  depends  only  on 
the  type  of  degrees  of  freedom  of  B. 

Example:  To  illustrate  these  concepts,  consider  the  case  of  an  arbitrary  (not 
necessarily  simply  connected)  ik-gon  B  translating  in  the  plane  amidst  a  collec- 
tion of  polygonal  obstacles  having  a  total  of  m  sides.  Each  placement  of  B  can 
be  specified  by  the  position  Z  of  some  fixed  reference  point  inside  B.  There 
are  n=0(km)  constraints  defining  FP,  and  they  induce  n  constraint  curves, 
where  each  such  curve  -y/  is  the  locus  of  placements  of  B  in  which  some 
specific  comer  of  B  touches  some  obstacle  edge,  or  some  specific  side  of  B 
touches  some  obstacle  comer.  In  this  special  case  it  is  easily  verified  that  each 
-y/  is  a  line  segment,  obtained  as  the  Minkowski  (vector)  difference  of  the 
obstacle  feature  and  the  feature  of  B  making  contact  (see  e.g.  [KLPS]  for 
more  details).  Thus  in  this  setting,  our  problem  is  to  calculate  the  face  in  an 
arrangement  of  n  line  segments,  which  contains  a  given  point  Zq  (see  [EGS], 
[EGHSSSW];  a  special  case  of  this  problem  has  also  been  studied  in  [PSS]). 
Such  a  face  is  illustrated  in  Figure  1. 

Let  r  be  the  collection  of  the  n  constraint  curves.  We  will  assume  that 
these  curves  are  in  general  position,  meaning  that  no  three  of  them  pass 
through  the  same  point,  that  no  two  are  tangent  to  one  another,  and  that  no 
two  of  them  overlap  (our  analysis  can  be  easily  modified  to  handle  degen- 
erate configurations  of  this  kind).  It  follows  that  any  pair  7/,  7^  of  curves  in 
r  intersect  in  at  most  d^  points  (By  Bezout's  Theorem).  (Actually  in  several 
cases  the  maximum  number  s  of  pairwise  intersections  can  be  shown  to  be 
considerably  smaller  than  that.)  Thus  the  total  number  of  intersections 
between  the  curves  in  F  is  at  most  sn  (n  - 1)/2  =  O(n^);  then  an  easy  applica- 
tion of  Euler's  formula  shows  that  the  arrangement  A(T)  of  the  arcs  in  T 
consists  of  C>(n^)  faces,  edges,  and  vertices. 
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Before  plunging  into  our  analysis,  let  us  digress  for  a  moment  to  con- 
sider several  special  cases  of  the  problem  in  which  the  complexity  of  the 
entire  space  FP  (not  just  of  a  single  component  of  it)  is  also  small.  To 
describe  these  cases,  suppose  for  a  moment  that  each  constraint  curve  -y<  is  a 
dosed  Jordan  curve,  and  that  it  partitions  the  plane  into  two  regions,  so  that 
one  of  them  (which  we  refer  to  as  the  interior  of  7/,  and  denote  it  by  Ki) 
consists  exclusively  of  forbidden  positions  of  fl  in  which  the  corresponding 
constraint  is  not  met,  whereas  in  the  complementary  region  that  constraint  is 
satisfied.  In  this  setting,  the  free  configuration  space  FP  is  just  the  comple- 
ment of  the  imion  AT  of  all  the  AT/'s.  If  the  maximum  number  of  intersections 
s  of  any  pair  of  curves  in  T  is  2,  then  it  was  shown  in  [KLPS]  that  the  overall 
combinatorial  complexity  of  K  (as  measured  e.g.  by  the  total  number  of 
intersections  of  the  "Vi's  which  lie  on  BK)  is  at  most  0{n).  A  very  recent 
result  [EGHPPSSS]  considers  the  case  in  which  each  Ki  is  the  region  enclosed 
between  the  ;c-axis  and  a  Jordan  arc  7/  whose  endpoints  lie  on  the  j:-axis,  so 
that  each  pair  of  the  arcs  yi  intersect  in  at  most  s  =  3  points,  and  shows  that 
in  this  case  the  overall  complexity  of  K  is  0(na(n)).  However,  as  soon  as 
J  >  4,  the  complexity  of  AT  can  become  ft(n^).  This  complexity  can  arise  in 
actual  2-DOF  motion  planning  problems,  e.g.  that  of  a  purely  translational 
motion  of  a  non-convex  polygonal  object,  or  that  of  a  horizontal  translational 
motion  of  a  vertical  line  segment  amidst  polyhedral  obstacles  in  3-space.  (It  is 
dear,  as  also  noted  above,  that  the  complexity  of  K  is  at  most 
0{sn^)  =  0(n2).) 

Thus,  unless  one  faces  particularly  favorable  spedal  cases,  the  complex- 
ity of  the  entire  FP,  and  thus  the  complexity  of  any  algorithm  that  computes 
the  entire  free  configuration  space,  must  be  quadratic  in  the  worst  case.  It  is 
the  purpose  of  this  p^er  to  show  that  the  combinatorial  complexity  of  a  sin- 
gle component  of  FP  is  in  general  close  to  linear,  and  that  calculation  of  such 
a  component  can  be  accomphshed  in  close  to  linear  time. 

Before  closing  this  initial  set  of  observations,  we  note  that  FP  can  be  cal- 
culated in  close  to  quadratic  time  in  the  worst  case: 

Theorem  2.1:  The  free  configuration  space  FP,  or,  more  generally,  the 
arrangement  A  of  the  n  constraint  curves,  can  be  calculated  in 
0{{n  +  p)  log  n)  =  0(n^log  n)  time,  where  p  is  the  number  of  intersections 
between  these  curves  yt. 

The  algorithm  which  calculates  FP  is  a  straightforward  modification  of  stan- 
dard line  sweeping  techniques  (see  [BO]),  appHed  to  the  collection  of  con- 
straint curves  7/.  (We  note  that  p  may  be  substantially  larger  than  the  actual 
number  q  of  these  intersections  which  he  on  dFP.  We  leave  it  as  an  open 
problem  whether  FP  can  be  calculated  in  time  0{{n+q)  log  n)  (or,  in  view 


of  the  results  below,  even  in  time  0{{n+q)  log^n))  where  ^  is  as  above.) 
Another  open  problem  is  whether  one  can  apply  the  topological  sweeping 
technique  of  [EG]  to  obtain  an  algorithm  whose  complexity  does  not  involve 
the  log  n  factor. 

Remark:  A  recent  result  of  [EGPPSS]  gives  an  0(nX,+2(n))  algorithm  for 
calculating  the  arrangement  A(r)  using  an  incremental  construction  tech- 
nique, whose  analysis  is  based  on  the  results  obtained  in  this  paper.  Another 
recent  result  of  [Q]  gives  a  randomized  algorithm  for  an  incremental  con- 
struction of  arrangements  of  line  segments,  whose  expected  time  complexity 
is  0(n  log  n  -I-  r),  where  t  is  the  number  of  intersections  between  the  given 
segments. 

3.  The  Complexity  of  a  Single  Component  of  FP 

In  this  section  we  obtain  an  almost-linear  upper  bound  on  the  combina- 
torial complexity  of  a  single  connected  component  of  FP.  Specifically,  we 
show 

Theorem  3.1:  Under  the  assumptions  made  in  the  preceding  section,  the  com- 
binatorial complexity  of  any  single  connected  component  of  FP  is  at  most 

Proof:  Let  /  be  the  given  connected  component,  and  let  C  be  a  connected 
component  of  its  boundary.  It  suffices  to  show  that,  if  k  arcs  of  F  appear 
along  C,  then  the  number  of  subarcs  of  these  arcs  which  constitute  C  is 
0{kj+2{k)).  Thus,  without  loss  of  generality,  we  may  assume  that  all  n  arcs 
of  r  appear  along  C.  For  each  7/  let  «/,  v,  be  its  endpoints.  Let  7/"*"  (resp. 
7/"  )  be  the  directed  arc  7^  oriented  from  ui  to  v^  (resp.  from  v/  to  w/). 

Without  loss  of  generality,  assume  C  is  the  exterior  boimdary  component 
of  /.  Traverse  C  in  coimterclockwise  direction  (so  that  /  lies  to  our  left)  and 
let  5  =  (ji,52.  •  •  •  >-s<)  be  the  circular  sequence  of  oriented  curves  in  F  in 
the  order  in  which  they  appear  along  C  (if  C  is  an  imbounded  arc,  5  is  a 
linear  rather  than  circular  sequence).  More  precisely,  if  during  our  traversal 
of  C  we  encounter  a  curve  7/  and  follow  it  in  the  direction  from  «/  to  v, 
(resp.  from  v/  to  «/)  then  we  add  7*  (resp.  7/")  to  5.  (As  an  example,  if  the 
endpoint  ui  of  7/  is  on  C  and  is  not  incident  to  any  other  arc,  then  traversing 
C  past  Ui  will  add  the  pair  of  elements  7/"  ,  7^  to  S,  and  symmetrically  for  v/. 
See  Figure  2  for  an  illustration.  Note  that  in  this  example  both  sides  of  an  arc 
7/  might  belong  to  our  connected  component;  in  our  original  motion  planning 
application  this  usually  was  not  the  case,  because  crossing  a  constraint  curve 
generally  means  that  the  system  is  passing  from  free  to  non-free  placements, 
but  the  generalized   problem  that  we  study  now   allows  for  such  two- 
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sidedness.) 

In  what  follows  we  will  use  the  followmg  notation.  We  denote  the 
oriented  arcs  of  T  as  ^i,  .  .  .  ,^2n-  For  each  ^,  we  denote  by  |4/|  the  non- 
oriented  arc  -y^  coinciding  with  ^/.  For  the  purpose  of  the  proof  we  will 
transform  each  arc  y,  into  a  very  thin  dosed  Jordan  curve  7*  by  taking  two 
non-intersecting  copies  of  7/  lying  very  dose  to  one  another,  and  by  joining 
them  at  their  endpoints.  This  will  perturb  the  component  / slightly  but,  imder 
OUT  assimiption  on  general  position,  it  will  not  change  the  combinatorial 
structure  of  the  boundary  of  /,  and  in  particular  of  C.  Note  that  this  transfor- 
mation allows  a  natural  identification  of  one  of  the  two  sides  of  7*  with  y^ 
and  the  other  side  with  7f  . 

We  next  need  the  following  lemmas: 

Lemma  3.2  (Ihe  Consistency  Lemma):  The  portions  of  each  arc  ^/  q^pear  in 
5  in  a  circular  order  which  is  consistent  with  their  order  along  the  oriented  ^/ 
(i.e.  there  exists  a  starting  point  in  S  (which  depends  on  ^/)  such  that  if  we 
read  S  in  circular  order  starting  from  that  point,  we  encounter  these  portions 
in  their  order  along  ^/). 

Proof:  Let  ^,  i)  he  two  portions  of  ^/  which  ^pear  consecutively  along  C  in 
this  order  (i.e.  no  other  portion  of  ^/  appears  along  C  between  ^  and  t]). 
Choose  two  points  x  i  I,  and  y  i  r\  and  connect  them  by  the  portion  a  of  C 
traversed  from  x  to  y,  and  by  another  arc  p  within  the  interior  of  \ii\* . 
Qearly  a  and  p  do  not  intersect  (except  at  their  endpoints)  and  they  are  both 
contained  in  the  complement  of  (the  interior  of)  /.  Thus  their  union  a  U  p  is 
a  closed  Jordan  curve  and  /  is  fully  contained  either  in  its  exterior  or  in  its 
interior.  We  claim  that  any  point  on  ^/  between  $  and  t]  is  contained  in  the 
side  of  aUP  which  does  not  contain  /.  Indeed,  connect  such  a  point  z  to  x 
along  an  arc  p  that  proceeds  very  near  4/  along  the  exterior  of  |^/ 1*  (see  Fig- 
ure 3).  Qearly  p  and  p  are  disjoint,  and,  deforming  a  slightly  as  necessary, 
we  can  assume  that  p  intersects  a  transversally  and  exactly  once,  which  is 
easily  seen  to  imply  our  claim.  This  daim  completes  the  proof  of  the  lemma. 

D 

For  eadi  directed  arc  ^,  consider  the  linear  sequence  V,  of  all  appearances 
of  4/  in  S,  arranged  in  the  order  they  appear  along  ^y.  Let  ji/  and  v,  denote 
respectively  the  index  in  5  of  the  first  and  last  element  of  V/.  Consider 
5  =  (ji,  .  .  .  ,J,)  as  a  linear,  rather  than  circular,  sequence  (this  step  is  not 
needed  if  C  is  unbounded).  For  each  arc  ^,,  if  jiy  >  v/  we  split  the  symbol  it 
into  two  distinct  symbols  ^,1,  ^/2,  and  replace  all  appearances  of  ^/  in  5 
between  the  places  \x.i  and  t  (resp.  between  1  and  v/)  by  ^,1  (resp.  ^/2).  (Note 
that  Lemma  3.2  implies  that  we  can  actually  spht  the  arc  ii  into  two  con- 
nected subarcs,  so  that  all  appearances  of  ^/i  in  5*  represent  portions  of  the 
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first  subarc,  whereas  all  appearances  of  4/2  represent  portions  of  the  second 
subarc.)  This  splitting  produces  a  sequence  5* ,  of  the  same  length  as  S,  com- 
posed of  at  most  4n  symbols. 

The  assertion  of  the  theorem  is  then  an  immediate  consequence  of  the 
following 

Lemma  3.3:  S*  is  an  (4n,s+2)  Davenport-Schinzel  sequence. 

Proof:  Since  it  is  dear  that  no  two  adjacent  elements  of  5*  can  be  equal,  it 
remains  to  show  that  5  *  does  not  contain  an  alternating  subsequence  of  the 
form  J  •  •  •  nq  ••  .  ^  •  •  •  -q  •  •  .  of  length  s  +  4.  Suppose  to  the  contrary  that 
S*  does  contain  such  an  alternation,  and  consider  any  four  consecutive  ele- 
ments of  this  alternation,  which,  without  loss  of  generality,  can  be  assumed 
to  be  C  •  •  •  r|  •  •  4  •  •  11-  Choose  points  x,y  €  i  <ind  points  z,w  €  t]  so  that 
C  passes  through  these  points  in  the  order  x,  z,  y,  w.  Consider  the  following 
five  Jordan  arcs: 

Pj^  =  an  arc  within  the  interior  of  |^  |*  connecting  x  to  y; 

P^  =  an  arc  within  the  interior  of  |ti  |*  connecting  z  to  w; 

P,j  =  the  portion  of  C  traversed  in  counterdodcwise  direction  from  a:  to  z; 

^jy  =  the  portion  of  C  traversed  in  counterclockwise  direction  from  z  toy; 

Pyn-  =  the  portion  of  C  traversed  in  counterdodcwise  direction  from  ytow. 

Note  that  p,,,  p^y,  p^  are  pairwise  non-intersecting  and  that  also  they 
do  not  intersect  pj^y,  p^^.  We  claim  that  p^^  and  p^^  must  intersect  one 
another.  Suppose  the  contrary,  and  consider  the  planar  graph  G  composed  of 
these  five  arcs  as  edges.  Qearly  G  has  three  faces,  and  its  edges  all  lie  in  the 
complement  of  the  interior  of  /.  Thus  /  is  fully  contained  in  just  one  of  these 
faces.  Moreover,  as  p  =  p,,  U  p^  U  p^^-  is  traced  from  j:  to  w,  all  points 
lying  on  the  left  side  of  p  suffidently  near  it  belong  to  /.  By  deforming  the 
plane  we  may  assume  that  p  lies  on  the  x  axis,  with  x,  z,y,w  appearing  along 
it  from  left  to  right  in  this  order,  that  the  portion  of  the  upper  half  plane  suf- 
fidently near  p  is  contained  in  /,  and  that  the  two  arcs  p^^  and  p^^.  emanate 
downwards  from  the  respective  points  y  and  z,  and  that  the  portions  of  these 
arcs  within  suffidently  small  neighborhoods  of  these  two  respective  points 
are  just  straight  vertical  segments;  see  Figure  4. 

Let  a  and  b  be  two  points  in  the  relative  interior  of  p^,,,  p^^.  respectively, 
let  a',  b'  be  two  corresponding  points  lying  below  and  very  close  to  a,  b 
respectively.  By  the  properties  of  G,  neither  a'  nor  b'  lies  in  the  face  of  G 
containing  /.  Moreover,  if  a'  and  b'  are  chosen  suffidently  near  p,  the 
straight  segment  connecting  a'  to  b'  will  intersect  each  of  ^^,  p^^  in  exactly 
one  point.  This  easily  implies  that  a'  and  b'  both  lie  in  the  same  face  <|)  of  G. 
Consider  the  closed  Jordan  curve  p*  traced  as  we  move  from  a  to  a'  along  a 
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short  segment,  from  a'  to  b'  along  a  path  that  connects  these  points  in  <J>, 
from  b'  to  b  along  another  short  segment,  and  finally  from  b  \o  a  along  I,.  It 
is  easily  checked  that  points  along  P,h.  slightly  after  z  lie  on  one  side  of  p* 
and  points  on  that  curve  just  before  w  lie  on  the  other  side  of  P* .  Thus  p^ 
has  to  intersect  p*,  which  however  is  impossible  because  p*  is  contained 
within  4>  U  p ,  which  is  openly  disjoint  from  p^^^  • 

This  shows  that  each  quadruple  of  consecutive  elements  in  our  alterna- 
tion induces  at  least  one  intersection  point  between  the  corresponding  arcs 
P;^  C  {  and  p^  C  Ti.  Moreover,  it  is  easily  diecked  that  for  any  pair  of  dis- 
tinct quadruples  of  this  type,  either  the  two  corresponding  subarcs  of  the 
form  P;jy  along  {  are  disjoint,  or  the  two  subarcs  p^^-  along  -r\  are  disjoint. 
Thus  all  these  intersections  must  be  distinct.  Since  the  number  of  such  qua- 
druples is  j  +  4  -  3  =  5  +  i,  we  obtain  a  contradiction,  which  completes  the 
proof  of  the  lemma,  and  thus  also  of  the  theorem,  n 

Remark:  It  has  been  recently  proved  in  [SS5]  that  if  the  curves  7/  in  F  are 
dosed  Jordan  curves,  or  Jordan  arcs  unbounded  in  both  directions,  then  the 
complexity  of  a  single  component  of  A(r)  is  at  most  0(X,(n)). 

Remark:  Applying  our  theorem  to  the  case  where  each  7/  is  a  line  segment, 
we  obtain  an  upper  bound  of  0{\-i{n))  =  0(na{n))  on  the  complexity  of  a 
single  component  of  A(r).  The  results  of  [WS]  and  [Sho]  imply  that  this 
bound  is  tight  in  the  worst  case.  This  special  case  of  the  theorem  has  already 
been  established  in  [PSS]. 

4.  Calculating  a  Single  Component  in  an  Arrangement  of  Curves 

In  this  section  we  show  how  to  calculate  a  single  component  /  in  an 
arrangement  of  a  collection  F  of  n  Jordan  arcs  71,  .  .  .  ,7„,  having  the  pro- 
perty that  no  pair  of  them  intersect  in  more  than  s  points.  Our  algorithm 
runs  in  time  0(X, +2(n)  log^n)  and  is  thus  close  to  linear.  We  assume  a 
model  of  computation  involving  infinite-precision  real  arithmetic,  in  which 
standard  operations  involving  one  or  two  curves  in  F  are  assumed  to  take 
constant  time.  Moreover,  since  our  technique  is  based  on  line  sweeping,  we 
assume  that  the  shape  of  each  curve  in  F  is  relatively  simple  and  not  too 
"wiggly".  Specifically,  we  assume  that  each  curve  in  F  has  at  most  t  points  of 
vertical  tangency,  for  some  fixed  constant  /,  so  that  we  can  break  it  into  at 
most  /  +  1  Jordan  arcs  that  are  monotone  in  the  x-direction.  Thus,  in  the 
remainder  of  this  section  we  will  assume  each  7  ^  F  to  be  j:-monotone  (this 
additional  condition  is  satisfied  in  most  applications;  in  particular,  it  holds  for 
curves  that  are  algebraic  of  a  fixed  maximal  degree).  Also,  we  assume  that 
the  curves  of  F  are  in  general  position,  so  that  each  intersection  of  a  pair  of 
these  curves  is  either  at  a  common  endpoint  or  is  a  transversal  intersection  at 
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a  point  in  the  relative  interior  of  both.  We  also  assiime  that  no  two  intersec- 
tion points  or  endpoints  lie  on  the  same  vertical  line,  so  as  to  simplify  the 
description  of  our  line  sweep  algorithm  (none  of  these  assumptions  are  essen- 
tial, and  simple  modifications  of  our  algorithm  will  make  it  work  also  in  the 
presence  of  degeneracies).  Typical  operations  that  are  assumed  to  take  con- 
stant time  are:  (a)  find  the  intersection  points  between  a  pair  of  curves  in  F; 
(b)  find  the  intersection  between  a  vertical  line  and  a  curve  in  F. 

The  algorithm  that  we  seek  thus  receives  as  input  a  collection  F  of  n  Jor- 
dan arcs  (or  dosed  curves)  with  the  above  properties,  and  a  point  x  not  lying 
on  any  of  these  curves.  Its  output  is  a  discrete  representation  of  the  con- 
nected component  /  of  the  arrangement  A  (F)  containing  x.  The  representa- 
tion that  we  will  use  is  the  collec^■on  of  the  connected  components  of  the 
boundary  of  /,  each  given  as  a  drcular  list  af  subaics  and  vertices  appearing 
along  that  boundary  component  in  counterclockwise  order.  The  algorithm 
extends  the  technique  in  [EGS]  for  the  calculation  of  a  component  (actually 
of  many  components  simultaneously)  in  an  arrangement  of  line  segments. 

The  high-level  description  of  our  algorithm  is  quite  simple.  We  use  the 
following  divide-and-conquer  technique.  We  spht  F  into  two  subcollections 
Fi,  F2  of  roughly  nil  curves  each,  calculate  recursively  the  components  /i, 
/2  of  A(Fi),  AiXi)  respectively,  that  contain  x,  and  then  "merge"  these  two 
components  to  obtain  the  desired  component  /.  Note  that  /  is  the  connected 
component  of /i  n  /2  containing  x.  However,  it  is  generally  too  expensive  to 
calculate  this  intersection  in  its  entirety,  and  then  select  the  component  con- 
taining X,  because  the  boundaries  of  /i  and  fi  might  intersect  in  many  (qua- 
dratically  many  in  the  worst  case)  points  that  do  not  belong  to  the  boundary 
of  /,  and  we  cannot  afford  to  find  all  of  them. 

The  set-up  for  the  merge  step  is  as  follows.  We  are  given  two  connected 
(but  not  necessarily  simply  connected)  regions  in  the  plane,  which  we  denote 
respectively  as  the  red  region  R  and  the  blue  region  B.  Both  regions  contain 
the  point  x  in  their  interior,  and  our  task  is  to  calculate  the  connected  com- 
ponent/of /?  n  B  which  contains  x.  The  boundaries  of  R  and  B  are  composed 
of  (maximal  connected)  portions  of  the  given  curves  in  F,  each  of  which  will 
be  denoted  in  what  follows  as  an  "arc"  (or  "subarc"). 

For  technical  reasons  that  will  be  explained  below,  we  extend  this  task  as 
follows.  Let  P  be  the  set  containing  x  and  all  endpoints  of  curves  of  F  which 
lie  on  the  boundary  of  either  Roi  B.  Qearly  P  contains  at  most  2n  + 1  points. 
For  each  w  €  P  let  /h,  denote  the  connected  component  of  ^  fl  fi  which  con- 
tains w  (these  components  are  not  necessarily  distinct).  Our  task  is  now  to 
calculate  all  these  components  (but  produce  each  distinct  component  just 
once,  even  if  it  contains  several  points  of  P).  We  will  refer  to  this  task  as  the 
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red-blue  merge.  (The  algorithm  given  below  actually  works  in  more  general- 
ity ~  R  and  B  can  be  the  union  of  several  connected  regions,  and  P  can  be  an 
arbitrary  (finite)  collection  of  points,  as  long  as  it  contains  all  the  endpoints 
of  the  curves  in  T  which  he  inside  B,  R  as  above.)  We  call  the  resulting  com- 
ponents /h,  purple  regions,  as  each  of  them  is  covered  by  both  the  red  and  the 
blue  region.  An  illustration  of  this  merge  is  shown  in  Figure  5. 

The  major  technical  result  on  which  our  algorithm  relies  is  that  the 
overall  complexity  of  these  purple  regions  is  small,  so  that  it  is  not  expensive 
to  produce  all  of  them.  This  will  be  a  consequence  of  the  following  extension 
of  the  combination  lemma  of  [EGS];  it  yields  a  somewhat  weaker  bound  than 
that  obtained  in  [EGS],  but  it  suffices  for  our  purpose. 

4.1.  The  combination  lemma  for  arrangements  of  curves 

The  combination  lemma  that  we  will  establish  in  this  subsection  is  some- 
what stronger  than  the  version  needed  for  our  red-blue  merge.  We  first 
introduce  a  few  notations.  Let  /?i,  .  .  .  ,/?^  be  a  collection  of  m  distinct  faces 
in  an  arrangement  of  a  collection  T^  of  "red"  Jordan  arcs,  and  let  fii,  .  .  .  ,5„ 
be  a  similar  collection  of  faces  in  an  arrangement  of  a  set  r^,  of  "blue"  Jordan 
arcs  (again,  each  pair  of  arcs  from  F,.  U  Fj,  are  assumed  to  intersect  in  at 
most  some  fixed  number  s  of  points).  Let  P  =  {pi,  .  .  .  ,pt}  be  a  collection 
of  points  so  that  each  Pi  i  P  belongs  to  one  red  face  R„^  and  to  one  blue  face 
5„,.  Let  El  be  the  connected  component  of  R„.  n  J9„.  containing  pi  (i.e.  £/  is 
the  face  of  the  combined  arrangement  of  F^.  U  Tf,  containing  p/).  Then  we 
have 

Lenuna  4.1  (The  Combination  Lemma  for  Arrangements  of  Curves):  The 
total  complexity  of  all  the  regions  £/  is  at  most  0(r+fc  +  ik),  where  r  and  b 
are  the  total  number  of  arcs  composing  the  boundaries  of  the  red  regions  and 
of  the  blue  regions  respectively. 

The  first  step  in  the  proof  is,  as  in  [EGS],  to  consider  the  special  case  of 
a  single  red  face  /?,  a  single  blue  face  B,  and  a  single  point  p  belonging  to 
B  r\R.  Let  E  be  the  connected  component  of  B  r\  R  containing  p.  Let  u 
(resp.  v)  denote  the  number  of  connected  components  of  BR  (resp.  BB),  and 
let  r  (resp.  b)  denote  total  complexity  (i.e.  number  of  subarcs  of  the  original 
curves  along  the  boundary)  of  R  (resp.  B).  We  assume  that  the  red  and  blue 
arrangements  are  in  general  position  so  as  to  avoid  degeneracies  in  the  struc- 
ture of  R,  B,  and  E,  such  as  tangencies  or  points  of  triple  intersection  of 
curves  along  their  boundaries.  This  assumption  is  made  for  the  sake  of  expo- 
sition; a  more  refined  version  of  the  analysis  given  below  can  handle  these 
degeneracies. 
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Lemma  4.2:  Under  the  non-degeneracy  assumptions,  the  complexity  of  £  is 
at  most  {s  +  3)(b  +  r  +  2u  +  2v  -At),  where  r  is  the  number  of  connected  com- 
ponents of  dE. 

Proof:  We  first  describe  the  general  outline  of  the  proof,  and  then  fill  in  the 
details  of  each  step.  We  analyze  each  component  4  of  the  boundary  of  E 
separately,  by  tracing  the  sequence  of  red  arcs  and  the  sequence  of  blue  arcs 
in  the  order  they  s^pear  along  t.  An  extension  of  the  "consistency  lemma" 
(Lemma  3.2)  shows  that  the  sequence  of  red  arcs  along  ^  is  consistent  with 
the  sequences  of  red  arcs  along  the  components  of  the  boundary  of  R,  and 
similarly  for  the  blue  sequence.  Notice  that  a  single  red  arc  from  R  can  be 
repeated  several  times  along  t,.  However,  these  repetitions  must  be  inter- 
spersed with  blue  arcs  which  "advance"  along  the  boundary  of  B.  Although  it 
is  possible  for  a  single  red  arc  to  be  interspersed  with  a  single  blue  arc  for  a 
while,  after  at  most  s+3  alternations  at  least  one  of  them  has  to  be  replaced 
by  another,  as  in  the  proof  of  Lemma  3.3.  Another  compUcation  that  can 
arise  is  that  I  may  visit  several  components  of  the  boundary  of  i?  or  of  B,  so 
that  some  of  them  are  visited  more  than  once.  We  show  that  the  duplication 
of  arcs  along  I  that  may  result  from  this  effect  is  only  linear  in  the  nxmiber  of 
components.  Altogether  these  argimients  imply  that  the  complexity  of  E  is 
linear  in  the  input  size,  as  asserted  in  the  lemma. 

To  begin  the  proof,  let  us  fix  a  single  component  ^  of  dE,  and  assume, 
without  loss  of  generality  that  it  is  the  exterior  component.  Trace  I,  in,  say, 
counterclockwise  direction  (with  E  lying  to  the  left),  and  let 
5  =  5{  =  (5i,  .  .  .  ,5„)  be  the  (circular)  sequence  of  the  subarcs  of  dR,  BB 
as  they  appear  along  ^.  Qearly  the  sum  of  the  lengths  of  5{,  over  all  con- 
nected portions  i  of  BE,  is  the  complexity  of  E. 

The  proof  consists  of  the  following  steps: 

(1)  Let  a  be  a  subarc  of,  say,  BR  which  appears  along  I,.  Let  a\,  ai  be 
two  connected  portions  of  aD^  consecutive  along  a,  such  that  when  a  is 
traversed  with  R  lying  to  its  left,  a^  precedes  ai.lx  follows  from  Lemma  3.2 
that  a  I  and  02  ^&  also  adjacent  along  {,  in  the  strong  sense  that  the  portion 
of  t  between  oi  and  ai  does  not  intersect  the  connected  component  of  BR 
containing  a. 

We  will  use  the  notation  "the  portion  of  S  between  si  and  s/'  to  mean  the 
subsequence  (5/  +  1,  .  .  .  ,jy_i)  if  i<j,  or  the  subsequence 
{si  +  i,  .  .  .  ,Sn,si,  .  .  .  ,sj-{)  if  j<i.  The  property  stated  above  then 
amounts  to  saying  that  if  Si=Sj  are  two  appearances  of  some  red  arc  a  in  5, 
with  the  first  appearance  preceding  the  second  one  along  a,  then  either  the 
portion  of  S  between  si  and  Sj  consists  of  blue  arcs  exclusively,  or,  if  it  con- 
tains red  arcs  at  all,  they  must  all  belong  to  other  components  of  BR. 
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(2)  Let  S^''^  be  the  (circular)  subsequence  of  5  obtained  by  deleting  from 
S  all  the  blue  subarcs,  and  let  S^''^  be  the  sequence  obtained  from  S^^^  by 
further  deleting,  from  left  to  right,  each  element  which  becomes  equal  to  the 
element  immediately  preceding  it.  The  sequences  S^''^  and  5^*^  are  defined 
symmetrically  for  the  blue  parts  of  5. 

We  claim  that  S^^^  is  of  length  at  most  r^  +  2u^-2,  where  «(  is  the 
number  of  distinct  connected  components  of  BR  appearing  along  ^,  and  r^  is 
the  total  number  of  red  subarcs  composing  these  «(  components. 

(Note  that  2  "{  ~  "»  ^°^  2  '"t  ~  '">  ^^^  do  component  of  BR  can 

appear  along  two  distinct  components  of  BE.) 

To  prove  the  claim,  suppose  a  is  a  red  subarc  appearing  more  than  once 
in  5'-''^  By  fl),  all  elements  of  S^''^  lying  between  two  consecutive  z^pear- 
ances  si  ,  sj  of  a  (arranged  in  this  order  along  a)  must  belong  to  other 
components  of  BR.  We  charge  the  second  appearance  of  a  to  the  component 
of  BR  containing  J/  +  i.  Let  ct^'')  be  the  (circular)  sequence  of  the  connected 
components  of  BR  in  the  order  they  appear  along  ^  (so  that  no  two  adjacent 
elements  of  ct^*")  are  equal).  As  in  the  proof  of  the  Q)mbination  Lemma  of 
[EGS],  it  can  be  shown  that  ct^*^)  is  a  circular  (u^,2)-  Davenport  Schinzel 
sequence  (i.e.  it  is  composed  of  uj  symbols,  no  two  adjacent  elements  of  it 
are  equal,  and  it  does  not  contain  a  subcycle  of  the  form 
(a  '  •  •  b  ■  •  •  a  •  •  ■  b)).  Hence  its  length  is  at  most  2«j— 2  (see  [ES]  for  a 
proof).  Moreover,  it  is  easily  checked  that  the  charging  scheme  described 
above  never  charges  an  element  of  o-^'')  more  than  once.  Hence  the  total 
number  of  duplications  of  elements  in  5^*^^  is  at  most  lu^—2,  from  which  the 
claim  follows. 

In  a  fully  symmetric  manner,  it  follows  that  the  length  of  5^*^  is  at  most 
fc{  +  2v{— 2,  where  vj,  b^  are  defined  analogously  as  the  number  of  con- 
nected components  of  BB  appearing  along  ^,  and  the  number  of  blue  subarcs 
composing  these  components. 

(3)  We  now  have  to  accoimt  for  duplications  of  adjacent  elements  in 
5^''^  5(*^  which  have  been  erased  in  5^'"\  5^*^  respectively.  Consider  S^^^ 
It  consists  of  (maximal)  rims,  where  each  run  is  a  contiguous  subsequence  of 
identical  elements,  and  is  represented  as  a  single  element  in  5^*^^  Let  p  be  a 
run  in  S^''^  of  length  /  of  the  single  subarc  a.  We  use  the  following  charging 
scheme.  For  each  element  jf*")  of  p,  other  than  the  first  element,  we  examine 
the  portion  8  of  5"^*)  which  appears  in  S  between  s^li  and  sY'> .  If  it  consists 
of  more  than  one  element,  then  they  must  all  be  distinct,  and  we  charge  one 
token  to  the  (element  of  S^^^  corresponding  to  the)  second  element  of  this 
portion.  Otherwise  8  consists  of  a  single  element  jj*^  If  sf"*  #  sf}\,  then 
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again  we  charge  one  token  to  the  cx)rrespondmg  element  of  5^*^ .  Otherwise 
no  charge  is  made.  We  use  a  symmetric  charging  scheme  to  the  runs  of  5^*^ . 

Consider  now  the  full  sequence  5.  An  element  of  S  is  said  to  be  charged 
if  it  is  a  charged  element  in  either  S^'^  or  5^*^ .  It  is  easy  to  check  that  no  ele- 
ment of  S^'^  or  of  S^*^  is  diarged  more  than  once,  so  the  total  number  of 
tokens  diarged  is  at  most 

Consider  a  portion  5  *  of  5  between  two  consecutive  diarged  elements  5/  and 
Sj.  Suppose  that  this  portion  contains  more  than  one  element  in  some  nm, 
and  suppose  that  the  rightmost  such  dupUcation  occurs  within  a  nm  of  some 
arc  a  in,  say,  S^''\  Let  s,^^=Sk^  =  a  be  this  rightmost  duphcation.  Since  no 
charge  was  made  within  5" ,  there  must  occur  a  single  blue  arc  x  between  5*^ 
and  Ji,,  which  furthermore  is  not  the  first  in  a  run  of  S^^K  Let  si^  =  si^=x  be 
the  two  consecutive  occurrences  of  x  within  its  run  in  5^*^  so  that  /2  =  ^i  +  l. 
If  j/j  also  Ues  within  S* ,  then  again  we  must  have  li  =  ki-l,  so  that  there  is 
a  single  appearance  of  a  between  these  two  appearances  of  x,  and  this  appear- 
ance is  also  not  the  first  in  its  run.  Continuing  backwards  in  this  manner,  it  is 
easily  verified  that  either  the  first  element  of  the  run  of  a  or  the  first  element 
of  the  nm  of  x  must  he  outside  (i.e.  before)  5* ,  and  that  a  and  x  are  the  only 
elements  being  dupUcated  in  consecutive  places  in  S^*"^  or  in  S^*^  within  S* . 
Now,  as  in  lemma  3.3,  the  maximum  number  of  alternations  of  the  arcs  a  and 
x  along  C  is  .y  -1-3,  so  that  at  most  s  +  2  alternations  can  occur  within  S* . 

In  other  words,  we  have  shown  that  the  excess  of  5^''^  and  5^*^  over  S^^^ 
and  5^*\  between  any  two  adjacent  charged  elements,  is  at  most  s  +  2.  Since 
the  number  of  charged  elements  is  at  most  r(  +  tj  +  2«j+2vj— 4,  it  follows 
that  the  total  length  of  5  is  at  most  {s  +  3)(r^  +  b^  +  2u^+2v^-4).  Summing 
over  all  components  i  of  dE,  we  obtain  that  its  total  complexity  is  at  most 
(j+3)(r  +  fe  +  2«+2v-4r),  as  asserted.  □ 

Remark:  Comparing  our  analysis  with  that  of  the  Combination  Lemma  of 
[EGS],  we  see  that  our  boimd  is  not  the  best  possible,  at  least  for  the  case  of 
straight  segments.  It  would  be  nice  to  sharpen  our  bound,  if  possible.  How- 
ever, our  result  imphes  that  the  complexity  of  E  is  0{r  +  b),  which  is  suffi- 
dent  for  our  purposes. 

Proof  of  the  combination  lemma:  The  proof  is  a  fairly  straightforward  (and 
somewhat  simplified)  ad^tation  of  the  proof  of  the  combination  lemma  for 
line  segments  given  in  [EGS].  For  the  sake  of  completeness  we  describe  the 
modified  proof  with  some  detail.  Fix  a  blue  face  B  =  Bj  which  contains  kj  of 
the  given  points,  say  pi,  .  .  .  ,pi,.  Let  Ci,  •  •  •  ,^/;  be  the  distinct  connected 
components  of  BB.  For  each  p,  let  £/  denote  the  connected  component  of 
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B  n  R,^  which  contains  p/,  as  defined  above.  Traverse  each  Cm  and  partition  it 
into  connected  portions  8  so  that  each  such  portion  intersects  the  boundary  of 
only  a  single  region  £/  (and  so  that  two  adjacent  portions  intersect  distinct 
such  regions);  note  that  in  general  the  (endpoints  of  the)  portions  8  are  not 
uniquely  defined.  We  define  a  plane  embedding  of  a  planar  graph  G  as  fol- 
lows. The  vertices  of  G  are  the  points  Pi,  .  .  .  ,Pkj  and  additional  Ij  points 
qi,  ...  ,qi  ,  so  that  q,  lies  inside  the  connected  component  Hi  of  R^  —  B 
whose  common  boundary  with  B  is  C/  •  For  eadi  portion  8  lying  on  some  l^ 
and  intersecting  some  BEi,  we  add  the  edge  {qm,Pi)  to  G,  and  draw  it  by  tak- 
ing an  arbitrary  point  in  8  fl  3£/  and  connect  it  to  p/  within  £/  and  to  q^ 
within  H„.  The  connectedness  of  each  £/  and  each  H„  implies,  as  in  [EGS], 
that  we  can  draw  all  edges  of  G  so  that  they  do  not  cross  one  another.  It  fol- 
lows from  the  definition  of  the  portions  6  that  in  this  embedding  of  G  each 
face  is  boimded  by  at  least  three  edges  (note  that  G  can  have  multiple  edges 
between  a  pair  of  vertices) .  Thus  by  Euler's  formula  the  number  of  edges  in 
G,  and  thus  the  number  of  portions  8,  is  at  most  3(kj  +  lj). 

We  next  define,  for  each  p/,  a  modified  "blue-red"  region  B*  containing 
Pi  as  follows.  U  R  =  R,^  does  not  intersect  dB  at  all,  then  we  take  B*  =  R. 
Otherwise  we  start  at  some  point  z  on  dB  D  BEi,  and  traverse  the 
corresponding  portion  8  of  dB  with  B  lying  to  the  left  imtil  its  last  intersec- 
tion with  dEt.  Then  we  turn  along  dR  into  B  and  follow  dR  imtil  its  next 
intersection  with  dB.  Since  this  intersection  necessarily  lies  in  dEi,  we  have 
landed  on  another  portion  8'  of  BB  which  intersects  BEt,  and  we  follow  8'  in 
the  direction  which  keeps  B  to  our  left,  until  its  last  intersection  with  BEi,  and 
continue  this  way  until  we  get  back  to  the  starting  point  z.  This  yields  one 
component  of  the  boundary  of  the  desired  blue-red  region  B* .  If  in  this  pro- 
cess we  have  not  encountered  all  portions  8  of  BB  intersecting  BEt,  we  pick 
another  starting  point  on  one  of  the  portions  we  have  missed,  and  repeat  the 
tracing  from  that  point.  Finally,  we  add  as  components  of  B*  all  components 
of  BR  which  boimd  £/  but  which  are  not  intersected  by  35  at  all.  Tracing  all 
components  of  the  boundary  of  BB*  in  this  way  yields  a  well  defined  con- 
nected region  bounded  between  these  boundaries. 

It  is  easily  checked  that  B*  contains  £,  and  is  contained  in  B.  We  define, 
in  a  completely  symmetric  manner,  a  modified  "red-blue"  region  R*  around 
each  pi.  It  follows  that  the  connected  component  of  B*  D  R*  which  contains 
Pi  is  exactly  Ef  Moreover,  assuming  that  no  two  points  pt,  Pj  give  rise  to  the 
same  intersection  face  £,  it  is  easy  to  check  that  BB*  D  Ej  =  BR*  H  Ej  =  0 
for  any  i¥=^j.  In  other  words,  no  arc  of  B*  or  of  R*  can  appear  on  the  boun- 
dary of  another  Ej  (with  Ej  lying  on  the  same  side  of  that  arc). 
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We  are  now  in  a  position  to  apply  Lemma  4.2.  Let 

bi  =  the  number  of  blue  arcs  of  B* , 

r,  =  the  number  of  red  arcs  of  R* , 

Ui  =  the  number  of  connected  components  of  dB* ,  and 

v/  =  the  number  of  connected  components  of  BR* . 

(By  construction,  all  such  components  actually  appeal  along  dEi.)  Since  each 
arc  of  dE(  is  either  a  (portion  of  a)  blue  arc  of  B*  or  a  (portion  of  a)  red  arc 
of  R* ,  a  slight  modification  of  the  proof  of  Lemma  4.2  (in  whidi  we  only 
need  to  accoimt  for  duplications  of  blue  arcs  of  B*  and  of  red  arcs  of  R*) 
implies  that  the  complexity  of  Et  is  at  most  0{bi  +  ri  +  u,  +  Vi).  Summing  these 
inequalities  over  all  points  pi,  we  conclude  that  the  overall  complexity  of  the 
regioi'S  £,  is  at  mvist 

0{^b,  +  ^n  +  2(«/  +  vO).  (i) 

/  /  / 

To  bound  ^  ^/.  consider  all  blue-red  regions  B*  contained  in  one  original 
blue  region  Bj .   Then    ^  bi  is  bounded  by  the  number  of  arcs  of  Bj  plus  a 

ptiBj 

term  proportional  to  the  number  of  subarcs  8  into  which  dBj  is  partitioned. 
By  the  preceding  argument,  this  additional  term  is  0{kj  +  lj)  where  kj  is  the 
number  of  points  pt  in  Bj,  and  Ij  is  the  number  of  connected  components  of 
dBj.  We  thus  obtain,  summing  over  all  blue  faces  Bj, 

^b,  =  0{b+  i  kj+  i  Ij). 

But  '^  kj  =  k,  and  the  total  number  of  components  of  the  blue  faces  is  clearly 

J 
bounded  by  their  total  complexity  b.  Hence 

j:b,  =  0{b  +  k).  (ii) 

Repeating  this  counting  for  the  red  faces  we  obtain 

^r,  =  0{r  +  k).  (iii) 

Finally,  ^  u/  (and  2  v/)  can  be  boimded  in  a  similar  manner,  noting  that  for 
each  original  blue  face  Bj ,  the  simi    ^    ui  is  bounded  by  Ij  plus  the  number 

PiiBj 

of  subarcs  8  along  dBj ,  so  that 

2  "/  =  0(  i  *y  +  i  'y )  =  Oib  +  k)  (iv) 

/  y=i  y=i 
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and  similarly 

2  V,  =  0(r  +  *)  .  (v) 

Combining  inequalities  (i)-(v)  completes  the  proof  of  the  lemma.  □ 

4.2.  The  red-blue  merge 

We  now  cx)ntinue  the  description  of  our  red-blue  merge.  Let  ir  be  the 
total  number  of  vertices  in  the  purple  regions.  The  Combination  Lemma  4.1 
implies  that  IT  =  0{b+r  +  n)  =  0(fc+r). 

To  facilitate  our  merge,  we  require  certain  information  to  be  precom- 
puted  and  available  for  each  collection.  Specifically,  we  require  that  the  red 
region  R  be  subdivided  into  j:-monotone  subregions  by  drawing  vertical  rays 
up  and  down  from  each  point  w  in  P  D  /?,  till  they  meet  an  edge  of  R  (we  call 
the  resulting  vertical  segment  through  w  the  red  vertical  divider  at  w,  and 
denote  it  by  pCw)).  It  is  easily  checked  that  this  does  produce  a  decomposi- 
tion of  /?  as  desired.  Similar  partitioning  is  required  for  the  blue  region, 
using  blue  vertical  dividers  (denoted  as  p(w)),  thereby  obtaining  a  similar 
collection  of  blue  j:-monotone  subregions. 

These  monotone  decompositions  of  the  red  and  blue  regions  are  easy  to 
obtain  using  a  straightforward  vertical  line  sweeping,  in  time 
0{{b  +  r)  log  (b  +  r)).  Note  that  a  particular  monotone  subregion  may  ter- 
minate on  the  left  or  the  right  either  because  of  a  point  of  P,  or  because  of  a 
locally  x-extremal  vertex  of  the  corresponding  region.  Our  algorithm  will 
produce  a  similar  partitioning  of  the  purple  regions  into  monotone  subre- 
gions, which  we  call  the  purple  subregions. 

We  calculate  the  purple  subregions  by  sweeping  with  a  straight  line. 
Notice  that  purple  subregions  start  or  end  at  j:-coordinates  associated  with 
either  a  point  of  P,  or  with  an  or-extremum  of  the  red  or  blue  region,  or  with 
a  red-blue  intersection.  In  a  left-to-right  sweep  we  will  discover  the  portion 
of  each  purple  subregion  that  is  to  the  right  of  the  leftmost  point  in  P  giving 
rise  to  it.  Then  afterwards,  in  a  right-to-left  sweep,  we  will  get  the  portion 
of  each  purple  subregion  to  the  left  of  the  rightmost  point  in  P  giving  rise  to 
it.  Together,  the  two  sweeps  discover  all  the  purple  subregions. 

Our  algorithm  will  thus  also  attempt  to  construct  purple  regions  incident 
to  each  blue  or  red  endpoint  (provided  that  the  endpoint  is  also  contained  in 
the  opposite-colored  region  -  otherwise  no  such  purple  region  is  to  be  gen- 
erated), even  though  such  a  purple  region  might  not  belong  to  the  desired 
purple  component  /  =  /,.  That  is,  some  of  these  purple  subregions  might  be 
disconnected  from  /  (see  Figure  5  for  an  example).  However,  we  will  be  able 
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to  detect  this  at  the  end  of  the  algorithm,  as  follows.  Consider  the  graph 
whose  nodes  are  the  final  purple  subregions,  and  whose  edges  connect  pairs 
of  purple  regions  adjacent  along  some  vertical  divider.  Then  a  simple 
breadth-first  search  on  this  graph,  starting  from  the  purple  subregion  that 
contains  x,  will  yield  the  desired  component  /;  all  the  unreachable  purple 
subregions  will  simply  be  discarded.  (Note  also  that  the  size  of  the  additional 
"fake"  purple  regions  is  at  most  0(^7  -I- r)  as  follows  easily  from  Lemma  4.1.) 

We  describe  below  only  the  left-to-right  sweeping  step;  the  right-to-left 
is  symmetric.  We  start  this  sweep  by  constructing  a  priority  queue,  ordered 
by  X  coordinate,  which  contains  all  the  vertices  appearing  along  the  boun- 
daries of  the  given  red  and  blue  regions,  together  with  the  point  x.  We  will 
be  sweeping  over  the  red  and  blue  subregions  separately,  and  at  the  same 
time  v'ill  also  he  sweepiug  over  ihe  purple  subregions,  ind  detect  (portions 
of)  them  as  we  sweep.  We  will  speak  of  the  red,  blue,  and  purple  planes 
respectively.  The  purpose  of  these  separate  sweeps  is  to  avoid  having  to  pro- 
cess "iminteresting"  red-blue  intersections,  i.e.  intersections  which  do  not 
occur  along  the  boimdary  of  a  purple  region. 

Every  time  we  encounter  a  point  p  of  P,  we  start  one  or  two  new  purple 
subregions  in  the  purple  plane.  At  such  an  event  we  create,  for  each  new  pur- 
ple region,  two  new  piirple  "scouts":  the  upper  scout,  and  the  lower  scout.  It 
is  the  job  of  these  two  scouts  to  walk  along  the  upper  and  lower  boundaries 
of  the  new  purple  regions,  respectively.  We  describe  the  behavior  of  the 
upper  scout  u  of  one  such  region;  the  lower  scouts  behave  symmetrically. 

The  upper  scout  u  starts  on  a  red  or  blue  arc,  as  determined  by  the 
dosest  of  the  upper  endpoints  of  the  vertical  dividers  p(p),  P(p).  The  scout  u 
will  move  right  along  that  arc  following  the  sweep  line,  but  it  needs  to  watch 
out  for  certain  events  that  might  influence  the  upper  purple  boundary  of  its 
region.  Without  loss  of  generahty,  we  assume  that  u  currently  sits  on  a  red 
arc. 

This  scout  u  has  first  of  all  to  look  up  to  the  nearest  blue  arc  ^  (note  that 
since  u  is  "purple",  the  entire  vertical  segment  between  u  and  p  must  be  con- 
tained in  B).  The  reason  is  that  the  blue  boundary  above  u  might  at  some 
future  point  drop  below  the  red  boundary  the  scout  is  currently  following.  If 
this  were  to  occur,  then  u  would  have  to  follow  the  blue  boundary,  because 
now  it  delimits  the  purple  subregion.  However,  there  might  be  another  scout 
V  already  watching  that  blue  arc  p  from  below.  In  that  case  only  the  highest 
of  u  and  v  needs  to  watch  P  from  below:  the  other  scout  can  rest,  since  it  is 
certain  that  the  higher  scout  is  "protecting"  it  from  p.  In  more  technical 
detail,  watching  a  blue  arc  p  means  that  u  has  to  determine  whether  the  arc  p 
it  lies  on  and  p  intersect  to  the  right  of  the  sweep  line,  and  if  so,  add  the 
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leftmost  such  intersection  as  an  event  into  the  priority  queue  of  the  sweep 
(each  of  these  operations  (except  for  the  priority  queue  insertion)  is  assumed 
to  take  constant  time  in  our  model  of  computation).  In  addition,  if  the  sweep- 
ing process  reaches  the  right  endpoint  of  either  p  or  p,  which  is  not  an  end- 
point  of  the  whole  curve  (from  F)  containing  that  arc,  u  has  to  retest  for  a 
potential  intersection  between  the  new  pair  of  arcs,  and,  if  it  exists,  add  the 
leftmost  such  intersection  to  the  queue.  When  such  a  red-blue  intersection  is 
eventually  swept  across,  u  has  to  move  to  the  blue  arc  p,  and  to  begin  to 
watch  the  red  boundary  above  it  (starting  with  p).  (Note  that  the  blue  boun- 
dary above  u,  or  the  red  boundary  it  currently  follows,  may  also  change 
discontinuously  if  the  sweep  reaches  the  endpoint  of  the  corresponding  full 
curve.  When  this  happens,  one  of  the  things  we  have  to  do  is  to  check 
whether  the  watching  assignment  of  some  purple  scout  has  to  change,  and  to 
inform  the  scout  of  this  change.  See  below  for  more  details.) 

The  scout  u  has  also  to  look  down  to  its  lower  partner  and  check  for 
their  (leftmost)  possible  intersection  (to  the  right  of  the  sweepline),  because 
when  the  two  of  them  come  together  the  current  purple  subregion  must  end. 
This,  of  course,  might  happen  earlier,  if  another  point  of  P  is  encountered 
between  these  scouts. 

The  key  property  here  is  that  each  blue  or  red  arc  is  watched  (at  any 
given  time)  by  at  most  one  upper  scout  and  at  most  one  lower  scout,  who  sit 
on  arcs  of  the  opposite  color.  But  these  assignments  of  who  watches  over 
whom  can  change. 

One  way  that  can  happen  is  that  a  purple  region  can  end  because  its  two 
purple  scouts  come  together.  This  will  occur,  for  instance,  when  the  right- 
most vertex  of  a  red  subregion  lies  inside  a  blue  subregion.  In  this  event,  the 
two  purple  scouts  of  that  purple  subregion  are  eliminated.  However,  some 
transfer  of  watching  responsibility  may  now  be  indicated.  If  the  upper  scout 
u  was  watching  a  blue  arc  p,  then  we  must  consult  the  next  upper  purple 
scout  down  from  u,  say  v.  If  v  is  currently  idle,  because  the  next  higher  blue 
arc  above  v  is  the  same  arc  p  watched  by  u,  then  u  transfers  to  v  the  responsi- 
bility of  watching  p.  If  v  is  already  watching  another  blue  arc,  then  we  leave 
it  imdisturbed,  as  its  blue  arc  must  be  below  p. 

Another  way  the  two  purple  partner  scouts  can  come  together  is  when  a 
purple  region  ends  at  a  red-blue  crossing.  Any  transfers  of  watching  responsi- 
bility that  need  to  happen  now  can  be  dealt  in  an  entirely  analogous  way.  If  a 
purple  region  ends  because  another  point  p  of  P  appears  between  the  scouts, 
then  again  the  two  scouts  are  eliminated,  but  in  this  case  they  will  generally 
be  replaced  by  new  scouts  spawned  by  p . 


23- 


The  reassignment  of  watching  responsibility  that  occurs  when  we  sweep 
through  a  point  p  of  P  is,  in  more  detail,  as  follows.  At  this  time  zero,  one  or 
two  new  purple  subregions  are  created.  Suppose  for  simphcity  that  only  one 
new  subregion  arises,  and  let  u  be  its  top  scout.  This  scout  u  finds  the 
opposite-colored  arc  e  it  has  to  watch  by  a  binary  search  through  the  hst  of 
arcs  of  that  color  currently  intersecting  the  sweepline.  But  then  u  also  has  to 
consult  the  upper  purple  scout  u  "^  (resp.  u "  )  lying  directly  above  u  (resp. 
below  m).  If  «"^  is  watching  the  same  e  (more  precisely,  if  u"^  lies  below  e) 
then  u  remains  idle.  Otherwise  u  begins  to  watch  e  and  checks  whether  u  ~  is 
also  watching  e,  in  which  case  u~  becomes  idle.  Similar  but  somewhat  modi- 
fied actions  are  taken  when  two  new  purple  regions  are  spawned  at  p . 

In  further  detail,  suppose  p  is  an  endpoint  of  a  red  original  curve;  then 
we  check  wneilicr  p  also  lie^  in  the  blue  region  B.  If  so,  we  start  one  or  t\^o 
new  purple  regions  incident  to  p  (and  lying  to  its  right),  and  proceed  with 
scout  creation  and  watching  reassignments  as  above;  otherwise  no  new  purple 
region  is  to  be  generated  at  p.  In  either  case,  if  the  vertical  divider  p(p) 
extends  both  up  and  down  from  p  (i.e.  p  is  a  point  of  "vertical  tangency"  on 
the  boundary  of  its  region),  then  some  red  boundaries  currently  watched  by  a 
purple  scout  may  change  discontinuously.  For  example,  if  the  two  red  arcs 
ei,  6  2  incident  to  p  extend  to  the  right  of  p  (with  ei  lying  below  €2),  then 
we  search  through  the  hst  of  purple  regions  along  the  sweepline  to  find  the 
upper  scout  u  lying  directly  below  p,  and,  if  it  lies  on  a  blue  arc,  tell  u  to 
start  watching  ei,  unless  it  is  watching  a  red  arc  lying  below  p.  If  «  lies  on  a 
red  arc,  then  if  it  is  idle  we  leave  it  imdisturbed;  if  it  is  watching  a  blue  arc 
lying  below  p  we  again  leave  it  undisturbed;  however,  if  u  watches  a  blue  arc 
P  above  p ,  then  it  is  easily  checked  that  p  must  also  lie  within  the  blue  region 
B,  and  consequently  two  new  purple  regions  will  be  created  at  p,  new  purple 
scouts  will  be  spawned  along  e-i,  62,  and  the  (upper)  scout  along  e^  will 
relieve  u  from  the  responsibihty  of  watching  p,  as  explained  above.  Similar 
or  symmetric  actions  are  taken  in  all  the  other  subcases. 

Finally,  some  transfer  of  watching  may  be  required  at  a  red-blue  crossing 
lying,  say,  on  the  top  boundary  of  some  purple  region.  Suppose  the 
corresponding  top  purple  scout  u  was  lying  on  a  red  arc  p  just  before  the 
intersection,  and  afterwards  it  moves  along  a  blue  arc  p.  As  noted  above,  u 
now  has  to  start  watching  p,  but  we  also  need  to  check  whether  the  top  pur- 
ple scout  V  lying  directly  below  u  has  to  change  the  arc  it  is  watching,  or 
become  idle.  Details  are  similar  to  the  cases  considered  above,  and  can  be 
easily  worked  out  by  the  reader. 

Note  that  this  procedure  not  only  produces  the  purple  regions,  but  also 
their  j:-monotone  decompositions  into  purple  subregions,  which  will  be  handy 
for  further  processing.  However,  some  preprocessing  might  still  be  required 
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before  subsequent  merges  (call  them  purple-violet)  can  be  performed.  This  is 
because  an  endpoint  p  of  some  purple  original  curve  may  also  lie  within  the 
violet  region,  in  which  case  we  will  need  to  find  the  violet  vertical  divider  at 
p,  and  this  information  in  general  is  neither  part  of  the  "purple"  data  nor  part 
of  the  "violet"  data,  and  can  be  obtained  only  by  combining  information  from 
these  two  collections  prior  to  their  merge. 

Let  us  now  analyze  the  complexity  of  this  process.  The  purple  scouts 
simply  trace  the  boundaries  of  the  purple  regions.  Each  such  scout  needs  to 
schedule  into  the  priority  queue  possible  intersection  events  between  the  arc 
it  is  currently  sitting  on,  and  the  arc  it  is  watching  (including  the  possible 
intersection  between  the  current  top  and  bottom  arcs  of  the  same  purple 
subregion).  Note  that  new  events  are  scheduled  when  we  sweep  either 
through  a  point  in  P,  through  a  blue  vertex,  through  a  red  vertex,  or  through 
a  red-blue  crossing  (which  is  a  vertex  of  a  purple  subregion).  Moreover,  at 
each  such  point  only  a  constant  number  of  new  events  are  scheduled.  Thus 
the  total  number  of  events  ever  scheduled  is  proportional  to  the  total  input 
and  output  size,  which,  by  Lemma  4.1,  is  0{b  +  r).  Thus  each  event  costs 
0(log  {b  +  r))  time  to  insert  into  (and  delete  from)  the  priority  queue.  The 
additional  operations  of  our  procedure  involve  updating  the  red,  blue,  and 
purple  lists  along  the  sweepline,  of  creating  and  eliminating  scouts  (i.e.  pur- 
ple subregions),  and  of  reassigning  watching  responsibihties.  It  is  plain  that 
we  need  to  perform  only  0{b  +  r)  such  operations,  and  that  each  can  be  car- 
ried out  in  0(log  {b  +r))  time  (because  the  maximimi  size  of  the  red,  blue, 
and  purple  lists  along  the  sweepline  is  at  most  0{b  +  r)).  The  final  step  of 
detecting  "true"  purple  subregions  (those  that  are  connected  to  the  "anchor" 
point  x),  and  of  eliminating  the  other  subregions,  can  be  done  by  a  simple 
gr^h  searching  (as  explained  earlier)  in  time  linear  in  the  number  of  purple 
subregions  produced  by  the  algorithm,  i.e.  in  0{b  +  r)  time.  Thus  our  pro- 
cedure runs  in  overall  time  0{{b  +  r)  log  (b  +  r)).  Moreover,  in  our  zqDplica- 
tion  each  of  the  regions  R  and  fl  is  a  single  connected  component  in  an 
arrangement  of  n/2  Jordan  arcs  (or  curves),  each  pair  of  whidi  intersect  in  at 
most  s  points.  We  have  thus  shown 

Theorem  4.3:  Given  two  connected  "red"  and  "blue"  regions  R,  B,  both  con- 
taining a  given  point  x,  whose  boundaries  are  composed  respectively  of  r  and 
b  subarcs  of  some  collection  of  Jordan  arcs  (or  curves),  no  pair  of  which 
intersecting  at  more  than  s  =  0{l)  points,  one  can  calculate  the  connected 
component  /  of  the  intersection  R  D  B  which  contains  x,  in  time 
0{{r  +  b)\og(r  +  b)),  under  the  assumptions  made  at  the  beginning  of  the 
section  concerning  the  given  curves  and  the  model  of  computation. 
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The  time  bounds  given  above  are  for  the  merge  step  of  our  algorithm. 
Using  Theorem  3.1,  we  thus  obtain,  by  straightforward  calculation, 

Theorem  4.4:  Given  a  collection  F  of  n  Jordan  arcs  (or  curves),  having  the 
property  that  no  pair  of  them  intersect  in  more  that  s  points,  and  also  satisfy- 
ing the  conditions  made  above,  one  can  calculate  the  connected  component  of 
A(T)  containing  a  specified  point  x,  in  time  C>(X,+2('»)  log^n)  (or  in  time 
0(X,(n)  log^n)  for  closed  Jordan  curves). 

Remarks:  (1)  This  result  follows  and  extends  the  previous  algorithm  given  in 
[EGS]  for  the  case  of  line  segments.  It  shows  that  the  red-blue  merge  is  a 
versatile  technique  of  a  purely  topological  nature,  which  can  be  applied  to 
fairly  general  classes  of  curves. 

(2)  Our  version  cf  the  Combinatiou  Lemma  (Lctnma  4.1)  is  weaker  than  the 
corresponding  lemmas  for  lines  or  for  line  segments,  as  given  in  [EGS],  in 
that  the  bound  it  produces  involves  the  red  and  blue  complexities  r,  b  with 
coefficients  greater  than  1  (namely  5-1-3).  This  does  not  allow  us  to  apply  the 
lemma  in  a  repeated  recursive  fashion  (as  has  been  done  in  [EGS])  to  obtain 
sharp  upper  boimds  on  the  complexity  of  many  connected  components  in  an 
arrangement  A(T)  as  above.  However,  assuming  that  a  sharp  bound  on  this 
complexity  can  be  obtained  by  other  means,  our  red-blue  merge  can  be  easily 
extended  to  obtain  an  algorithm  for  calculating  m  such  components,  each 
specified  by  a  given  point  within  it,  in  time  comparable  with  their  total 
worst-case  complexity.  Recently,  [CEGSW]  have  obtained  sharp  upper 
bounds  for  the  complexity  of  m  faces  in  arrangements  of  n  circles,  or  of  n 
imit  circles,  or  of  n  pseudo-lines  (i.e.  a:-monotone  imbounded  arcs,  each  pair 
of  which  intersects  at  most  once).  Using  their  bounds,  we  can  obtain  the  fol- 
lowing results  (we  omit  details  of  the  solutions  of  the  resulting  divide-and- 
conquer  recurrences,  which  are  based  on  random  sampling  of  the  given  col- 
lection of  curves,  and  are  very  similar  to  those  obtained  in  [EGS]). 

Corollary  4.5:  (a)  One  can  calculate  m  distinct  faces  in  an  arrangement  of  n 
unit  circles  in  (randomized)  time  0(m^'^"*n^'^"*"^'  log  n  -I-  n  log^/i),  for  any 
8>0. 

(b)  One  can  calculate  m  distinct  faces  in  an  arrangement  of  n  arbitrary  circles 
in  (randomized)  time  0{m^'^~^n*'^'^^^  log  n  -I-  n  log^n),  for  any  8>0. 

(c)  One  can  calculate  m  distinct  faces  in  an  arrangement  of  n  pseudolines 
(under  an  ^propriate  model  of  computation)  in  (randomized)  time 
0(m2/3-8„2/3+28  log  „  +  „  log^n),  for  any  8>0. 

(3)  Concerning  lower  bounds,  it  is  easy  to  establish  an  Cl{n  log  n)  bound  by 
reducing  from  sorting.  We  do  not  know  of  any  larger  lower  bound.  In  partic- 
ular, we  pose  it  as  an  open  problem  whether  a  component  of  A(r)  can  be 
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computed  in  time  0(X,+2('»)  log  ")•  We  note  that  an  algorithm  with  such 
complexity  exists  for  the  calculation  of  the  lower  envelope  of  the  arcs  in  T 
([At],  [HS]). 

(4)  If  the  collection  T  consists  of  n  closed  Jordan  curves,  then  we  can  use  our 
red-blue  merge  to  calculate  a  single  component  of  A{T)  in  time 
0(kf{n)  log^n),  making  use  of  the  bound  obtained  in  [SS5],  as  mentioned  at 
the  end  of  section  3.  In  the  special  case  5 =2,  a  simple  adaptation  of  the  algo- 
rithm in  [KLPS]  yields  an  algorithm  with  that  complexity. 

(5)  Returning  to  the  original  motion  planning  problem,  the  above  procedure 
will  yield  a  single  connected  component  of  FP,  provided  this  component  is 
fully  contained  in  a  single  planar  patch  of  the  entire  parametric  space  AP.  If 
AP  is  not  planar,  and  the  desired  component  C  which  contains  the  initial 
placement  Z  "spills"  over  into  more  than  one  patdi,  we  obtain  C  using  the 
following  technique.  Let  APi ,  .  .  .  ,AP,  be  the  planar  patches  which  compose 
AP,  and  suppose  that  Z  lies  in  APi.  We  first  calculate  the  component  Co  in 
APi  containing  Z.  If  Co  does  not  reach  the  boundary  of  APi  then  C  =  Co- 
Otherwise,  we  also  calculate  all  the  unboimded  components  in  each  of  the 
patches  APj.  Assuming  that  the  number  of  patches  is  constant,  and  that  each 
of  the  constraint  curves  is  broken  into  a  constant  number  of  connected 
subarcs,  each  lying  within  a  single  patch,  it  is  easUy  checked  that  the  total 
complexity  of  all  unbounded  components  of  the  patches  is  0(X,+2(«))>  and 
that  they  can  all  be  calculated  in  time  0{Kj+2in)  log^n),  using  the  above 
algorithm.  The  desired  component  C  is  easily  seen  to  be  the  union  of  Co 
with  certain  portions  of  the  unbounded  components  within  the  patches. 
More  precisely,  as  Co  reaches  the  boimdary  of  APi,  C  extends  into  another 
patdi,  necessarily  as  a  portion  of  the  unbounded  component  within  that 
patch;  similar  "propagation"  of  C  into  further  patches  continues  imtil  the 
whole  of  C  is  obtained.  Qearly  this  tracing  of  C  can  be  accompUshed  within 
the  above  time  bound. 

Let  us  discuss  a  few  appUcations  of  our  results  to  some  specific  motion 
planning  problems.  As  a  first  example,  consider  a  line  segment  B  =  PQ  free 
to  rotate  in  three  dimensions  about  its  endpoint  P  which  is  fixed  at  the  origin. 
The  motion  of  B  clearly  has  two  degrees  of  freedom,  and,  assuming  the  obs- 
tacles that  B  has  to  avoid  are  all  polyhedral  (with  a  total  of  n  edges),  one 
easily  diecks  that  any  pair  of  constraint  curves  intersect  in  at  most  j  =  1  point. 
Thus  the  complexity  of  a  single  component  of  FP  is  0{na{n)),  and  can  be 
calculated  in  time  0(na(n)  log^n).  (Note  that  in  this  case  AP  can  be 
represented  as  the  surface  of  a  sphere  in  three  dimensions,  or,  alternatively, 
as  the  union  of  two  planar  patches.) 
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As  a  second  example,  let  B  be  a  two-link  arm  PQR  moving  in  the  plane 
(amidst  polygonal  obstacles)  with  the  point  P  fixed  at  the  origin.  Here  AP  can 
be  represented  by  a  torus,  or  by  a  union  of  four  planar  patches.  By  using 

tan  —,  tan  -;^  as  the  system  parameters,  it  is  easily  checked  that  each  con- 

straint  curve  is  algebraic  of  (at  most)  fourth  degree.  Moreover,  it  can  be 
shown  that  each  pair  of  constraint  curves  intersect  in  at  most  s=2  points  (on 
a  single  patch).  It  follows  that  the  complexity  of  a  single  component  of  FP  is 
0(X4(n)),  and  that  it  can  be  calculated  in  time  0(X4(n)  log^n). 

As  a  final  example,  consider  the  problem  of  coordinated  motion  planning 
for  two  "planar  Stanford  arms"  [FWY],  [SSi],  where  each  arm  is  a  straight 
segment  PiQ,  free  to  translate  through  a  fixed  point  C,  and  also  rotate  about 
that  point,  for  i  =  l,2.  We  consider  a  spedai  case  of  coordinated  motion,  in 
which  the  two  endpoints  d,  Q2  must  be  in  contact  with  each  other 
throughout  the  motion.  This  problem  has  been  studied  in  [FWY],  who  gave 
an  0(n^  log  n)  algorithm  (where  n  is  the  total  number  of  obstacle  comers  and 
edges).  A  recent  result  in  [SSi]  shows  that  the  desired  free  space  can  be  taken 
to  be  a  single  component  in  the  intersection  of  the  two  2-D  configuration 
spaces  of  the  two  separate  arms.  Since  each  of  these  spaces  has  0(n)  com- 
plexity [FWY],  [SSi],  a  single  appUcation  of  our  red-blue  merge  yields  the 
desired  space  in  time  0{n  log  n). 

S.  Conclusion 

The  results  obtained  in  this  paper  provide  a  satisfactory  solution  to  the 
general  motion  planning  problem  with  two  degrees  of  freedom.  The  problem 
is  in  a  much  more  confused  state,  however,  when  we  consider  problems  with 
three  degrees  of  freedom.  There,  under  assumptions  similar  to  those  made 
above,  the  complexity  of  the  entire  FP  is  O(n^)  in  general,  although  there 
are  certain  favorable  cases  in  which  the  complexity  reduces  to  quadratic  or 
near-quadratic.  For  example,  for  a  line  segment  moving  in  the  plane  amidst 
polygonal  obstacles,  the  complexity  of  FP  is  0{n^)  [LS2],  [SiS].  For  a  convex 
k'gon  moving  in  a  similar  environment,  the  complexity  of  FP  is  0{kn\s{kn)) 
[KS2].  However,  for  a  non-convex  (e.g.  an  L-shaped)  polygon,  the  complex- 
ity of  FP  can  be  ft(n^),  as  is  easily  checked.  As  above,  we  would  like  to  con- 
jecture that  a  single  connected  component  of  FP,  which  is  all  we  really  need 
to  compute,  will  have  smaller,  e.g.  near-quadratic,  complexity.  However  this 
appears  to  be  a  much  harder  problem,  and  is  largely  open.  Some  progress 
was  recently  made  on  this  problem  for  the  special  case  in  which  the  surface 
patches  bounding  FP  are  n  triangles  in  3-space.  It  was  shown  in  [PS]  that  the 

3-  — 
complexity  of  a  single  component  of  their  complement  is  at  most  0(n     ^'  ). 
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This  was  improved  in  [AS]  to  0(n'"^*')  for  any  €>0,  where  this  also  bounds 
the  total  complexity  of  all  non-convex  components.  Another  simplified  case 
is  where  instead  of  calculating  the  boundary  of  a  component  of  the  comple- 
ment of  the  given  constraint  surfaces,  we  only  want  to  calculate  their  upper 
or  lower  envelope,  i.e.,  the  portions  of  these  surfaces  seen  from  a  point  at 
infinity  in  the  direction  of  the  positive  or  negative  z-axis.  For  triangles,  it  was 
shown  in  [PS]  that  the  complexity  of  their  upper  envelope  isO{n^a(n)),  and 
that  this  bound  is  worst-case  optimal.  However,  for  arbitrary  surfaces,  even 
in  this  restricted  problem  subcubic  upper  bounds  are  not  known. 
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Figure  1.  A  face  in  an  arrangement  of  segments. 
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Figure  2.  Traversing  a  face  boundary  and  the  resulting  sequence  S. 
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Figure  3.  Dlustration  of  the  proof  of  the  consistency  lemma. 
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Figure  4.  Dlustration  of  the  proof  of  Lemma  3.3. 


Figure  5.  The  red-blue  merge. 
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